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Bose-Einstein condensation in a one-dimensional interacting system due to power-law 

trapping potentials 
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We examine the possibility of Bose-Einstein condensation in one-dimensional interacting Bose gas 
subjected to confining potentials of the form V ex _t(x) = Vo(|a;|/a) 7 , in which 7 < 2, by solving 
the Gross-Pitaevskii equation within the semi-classical two-fluid model. The condensate fraction, 
chemical potential, ground state energy, and specific heat of the system are calculated for various 
values of interaction strengths. Our results show that a significant fraction of the particles is in the 
lowest energy state for finite number of particles at low temperature indicating a phase transition 
for weakly interacting systems. 
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I. INTRODUCTION 

The recent observations of Bose-Einstein condensation 
(BEC) in trapped atomic gases |l]-|5| have renewed inter- 
est in bosonic systems H0. BEC is characterized by a 
macroscopic occupation of the ground state for T < To, 
where To depends on the system parameters. The success 
of experimental manipulation of externally applied trap 
potentials bring about the possibility of examining two or 
even one-dimensional Bose-Einstein condensates. Since 
the transition temperature To increases with decreasing 
system dimension, it was suggested that BEC may be 
achieved more favorably in low-dimensional systems ||. 
The possibility of BEC in one - (ID) and two-dimensional 
(2D) homogeneous Bose gases is ruled out by the Hohen- 
berg theorem ||. However, due to spatially varying po- 
tentials which break the translational invariance, BEC 
can occur in low-dimensional inhomogeneous systems. 
The existence of BEC is shown in a ID noninteracting 
Bose gas in the presence of a gravitational field Jl(J, an 
attractive-^ impurity pT| , and power-law trapping po- 
tentials |l2j . Recently, many authors have discussed the 
possibility of BEC in ID trapped Bose gases relevant 
to the magnetically trapped ultracold alkali-metal atoms 
|l3| |l8| . Pearson and his co-workers (ll]] studied the in- 
teracting Bose gas in ID power-law potentials employing 
the path-integral Monte Carlo (P1MC) method. They 
have found that a macroscopically large number of atoms 
occupy the lowest single-particle state in a finite system 
of hard-core bosons at some critical temperature, ft is 
important to note that the recent BEC experiments are 
carried out with finite number of atoms (ranging from 
several thousands to several millions), therefore the ther- 
modynamic limit argument in some theoretical studies 
[jrj| does not apply here || . 

The aim of this paper is to study the two-body interac- 
tion effects on the BEC in ID systems under power-law 
trap potentials. For ideal bosons in harmonic oscillator 
traps transition to a condensed state is prohibited, ft 
is anticipated that the external potentials more confin- 



ing than the harmonic oscillator ty pe would be possible 
experimentally, ft was also argued jl5j that in the ther- 
modynamic limit there can be no BEC phase transition 
for nonideal bosons in ID. Since the realistic systems are 
weakly interacting and contain finite number of particles, 
we employ the mean-field theory (2(]j2l|] as applied to a 
two-fluid model. Such an approach has been shown to 
capture the essential physics in 3D systems [ pi] ]. The 
2D version [£2| is also in qualitative agreement with the 
results of PIMC simulations on hard-core bosons p3[ . 
In the remaining sections we outline the two-fluid model 
and present our results for an interacting ID Bose gas in 
power-law potentials. 



II. THEORY 

In this paper we shall investigate the Bose-Einstein 
condensation phenomenon for ID interacting Bose gas 
confined in a power-law potential: 



V ext (x) = V ^ 
1 a 



(1) 



where Vo and a are some suitable energy and length pa- 
rameters defining the external potential, and 7 controls 
the confinement strength. Presumably, they can be ex- 
perimentally adjusted. Using the semi-classical density of 
states, the transition temperature To and the fraction of 
condensed particles N /N for the noninteracting system 
were calculated as fpJ 
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where k = 2(2m) ' a/jh V {in is the mass of bosons 
and h is the Planck's constant), and 



1 



T( 7 ) 



1 x 1 ^- 1 dx 



(4) 



and 



roc l/ 7 -l/2 i 

G(7) = / ^j— = r(l/ 7 + 1/2) C(l/7 + 1/2) , 

Jo e — i 

(5) 

in which T(x) and £(x) are the gamma and the Riemann 
zeta-functions, respectively. The total energy of the non- 
interacting system for T <Tq (/j = 0) is given by 



r(l/ 7 + 3/2) C(l/7 + 3/2) / T \ 1/7+3/2 



(E) 

Nk B T r(l/ 7 + 1/2) C(l/7 + 1/2) \T 



(6) 



Figure 1 shows the variation of the critical temperature 
To as a function of the exponent 7 in the trapping po- 
tential. It should be noted that To vanishes for harmonic 
potential due to the divergence of the function G(j = 2). 
It appears that the maximum To is attained for 7 0.5, 
and for a constant trap potential (i.e. V ex t(%) = Vo) the 
BEC disappears consistent with the Hohenberg theorem. 
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FIG. 1. The variation of the critical temperature To with 
the external potential exponent 7. 

We are interested in how the short-range interaction 
effects modify the picture presented above. To this end, 
we employ the mean-field formalism and describe the col- 
lective dynamics of a Bose condensate by its macroscopic 
time-dependent wave function T(x, t) — ^(x) exp (—i/it), 
where /i is the chemical potential. The condensate wave 
function ty(x) satisfies the Gross-Pitaevskii (GP) equa- 
tion MM 



fir d 2 

- — -j— «■ + V ext (x) + 2gm(x) + gV 2 (x) 
2m dx z 



where g is the repulsive, short-range interaction strength, 
and ni (x) is the average noncondensed particle distribu- 
tion function. We treat the interaction strength g as a 
phenomenological parameter without going into the de- 
tails of actually relating it to any microscopic descrip- 
tion [P6[ . In the semi-classical two- fluid model |27 28 
the noncondensed particles can be treated as bosons in 
an effective potential |2l|,|29|| 



V cS (x) = V ext (x) + 2gm{x) + 2g^ 2 (x) 
The density distribution function is given by 
, v f dp 1 



(8) 



2Trh exp{[p 2 /2m + V eff (x) - fi]/k B T} - 1 ' 

(9) 

and the total number of particles TV fixes the chemical 
potential through the relation 



N = N + 



p{E) dE 



exp [(£7 - ix)/k B T] - 1 



(10) 



where No = J V P 2 (x) dx is the number of condensed parti- 
cles, and the semi-classical density of states is determined 
by 



p(E) 



dx 



V B a(x)<E \JE - V c s(x) 



(11) 



The GP equation yields a simple solution when the ki- 
netic energy term is neglected (the Thomas-Fermi ap- 
proximation) 



* 2 0e) 



M ~ Vcxtjx) - Zgnijx) 
9 



% - V cxt (x) - 2gn 1 (x)) , 
(12) 



&(x) = fj$(x) , 
(7) 



where 0[x] is the unit step function. More precisely, the 
Thomas-Fermi approximation J7 pc|j30| ] would be valid 
when the interaction energy ~ gNo/K, far exceeds the 
kinetic energy fi 2 /2mA 2 , where A is the spatial extent 
of the condensate cloud. For a linear trap potential 
(i.e. 7 = 1), a variational estimate for A is given 

by A = [fi 2 /2m(TT/2) 1 ^2a/V ] 1/3 . We note that the 
Thomas- Fermi approximation would breakdown for tem- 
peratures close to To where No is expected to become 
very small. 

The above set of equations [Eqs. (9)-(12)] need to be 
solved self-consistently to obtain the various physical 
quantities such as the chemical potential fi(N,T), the 
condensate fraction Nq/N, and the effective potential 
V e f[. In a 3D system, Minguzzi et al. [Jn) solved a simi- 
lar system of equations numerically and also introduced 
an approximate semi-analytical solution by treating the 
interaction effects perturbatively. Motivated by the suc- 
cess [plp2]| of the perturbative approach we consider a 



2 



weakly interacting system in ID. To zero-order in gn\(r), 
the effective potential becomes 

if fi < U cxt (x) 



N = N Q + K(k B T) 1 ^+ 1 / 2 J(7,M,r) , (16) 



where 

v ( „\ _ / V ext (x) if fJ,<V cxt (x) n „, 

1/cffW -\2 M -Kxt(x) if M>Kxt(x). ^ /( 7)/U ,T) = F( 7 ) / 



x 1 /"*- 1 / 2 dx 
ze x — 1 



Figure 2 displays the typical form of the effective po- 
tential within our semi-analytic approximation scheme. 
The most noteworthy aspect is that the effective poten- 
tial as seen by the bosons acquire a double-well shape 
because of the interactions. We can explain this result 
by a simple argument. Let the number of particles in 
the left and right wells be Nl and Nr, respectively, so 
that N = Nl + Nr. The nonlinear or interaction term 
in the GP equation may be approximately regarded as 
V = N 2 + N R . Therefore, the problem reduces to the 
minimization of the interaction potential V, which is 
achieved for Nl = Nr. 
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FIG. 2. Effective potential V e g(x) in the presence of in- 
teraction (xq = (/i/Fo) 1 ^ 7 cs). Thick dotted line represents 
external potential F ex t(a:)- 

The number of condensed atoms is calculated to be 
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The density of states is given by 



(14) 



p(E) = k 



H(j,n,E) (2/j, - E) 1 /^ 1 / 2 if fi<E<2[i 



F( 7 ) £V7-i/2 



if E>2p, 



(15) 



where 



E/^-E) x l h -l dx 



Using the above density of states, conservation of total 
number of particles gives us a transcendental equation 
for the chemical potential 



2ft/k B T 

,^ /kBT TT , , ^ (2^/k B T-x) 1 ^- 1 / 2 dx 



ze x — 1 



in which z = e~v-l kBT . The chemical potential /i(N,T) 
is determined from the solution of Eq. (16). Finally, the 
total energy of the interacting system can be written as 

(E) = [(E) nc (N-N )/2 + (E) c ]/N , (17) 

where (E) nc is the energy of the noncondensed particles 

f Ep(E) dE 

[ Uc J cxp [(E — p)/k B T] — 1 

= K{k B Tfh+V2 j( 7 , M ,T), (18) 



where 



J2 



x l/ 7 +l/2 dx 

ze x — 1 



2fi/k B T 

2 »' k * T , (2 t i/k B T-x) 1 l"'+ 1 l 2 dx 
(7 ' ^ X) z~^~l • 

and (E) c is the energy of the particles in the condensate 

2a 7 V +1/7 



^)c = f/ dx = - 



(i + 7 )(2 7 + i) ff u J 



l/ 7 



• (19) 



The kinetic energy of the condensed particles is neglected 
within our Thomas- Fermi approximation to the GP equa- 
tion. 



III. RESULTS AND DISCUSSION 

Up to now we have based our formulation for arbitrary 
7 , but in the rest of this work we shall present our re- 
sults for 7 = 1. Our calculations show that the results 
for other values of 7 are qualitatively similar. In Figs. 
3 and 4 we calculate the condensate fraction as a func- 
tion of temperature for various values of the interaction 
strength r\ = g/Voa (at constant N = 10 5 ) and different 
number of particles (at constant r) = 0.001), respectively. 
We observe that as the interaction strength r\ is increased, 
the depletion of the condensate becomes more apprecia- 
ble (Fig. 3). As shown in the corresponding figures, a 
significant fraction of the particles occupies the ground 
state of the system for T < Tq. The temperature depen- 
dence of the chemical potential is plotted in Figs. 5 and 
6 for various interaction strengths (constant N — 10 5 ) 
and different number of particles (constant r? = 0.001) 
respectively. 
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FIG. 5. The temperature dependence of the chemical 
FIG. 3. The condensate fraction No/N versus temperature potential fJ,(N,T) for various interaction strength and for 
T/To for N = 10 5 and for various interaction strengths 77. N = 10 5 particles. 



Effects of interactions on /z(iV, T) are seen as large de- 
viations from the noninteracting behavior for T < Tq. In 
Fig. 7 we show the ground state energy of an interacting 
ID system of bosons as a function of temperature for dif- 
ferent interaction strengths. For small 77, and T < Tq, 
(E) is similar to that in a noninteracting system. As 
77 increases, some differences start to become noticeable, 
and for 77 sa 1 we observe a small bump developing in 
(E). This may indicate the breakdown of our approxi- 
mate scheme for large enough interaction strengths, as we 
can find no fundamental reason for such behavior. It is 
also possible that the Thomas-Fermi approximation em- 
ployed is violated as the transition to a condensed state 
is approached. 




Although it is conceivable to imagine the full solution 
of the mean-field equations [Eq. (9)-(12)] may remedy the 
situation for larger values of 77, the PIMC simulations 
also seem to indicate that the condensation is inhibited 
for strongly interacting systems. The results for the spe- 
cific heat calculated from the total energy curves, i.e. 
CV = d(E)/dT, are depicted in Fig. 8. The sharp peak 
at T — T tends to be smoothed out with increasing in- 
teraction strength. It is known that the effects of finite 
number of particles are also responsible for such a be- 
havior p0[ |. In our treatment these two effects are not 
disentangled. It was pointed out by Ingold and Lam- 
brecht |l4| that the identification of the BEC should also 
be based on the behavior of CV around T ~ T . 
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50 ■ 




-50 1 ■ 1 ' 1 ■ 1 ' ' 1 1 1 1 

0.0 0.2 0.4 0.6 0.8 1.0 1.2 



FIG. 4. The condensed fraction N /N versus temperature FIG - 6 ' The temperature dependence of the chemical po- 

T/T for 77 = 0.001 and for different number of particles N. tential P( N < T ) for different number of Particles N and for 

77 = 0.001. 
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FIG. 7. The temperature dependence of the total energy of 
ID Bose gas for various interaction strengths r\ and TV = 10 5 
particles. 

Our calculations indicate that the peak structure of CV 
remains even in the presence of weak interactions, thus 
we are led to conclude that a true transition to a Bose- 
Einstein condensed state is predicted within the present 
approach. 
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bosons. Thus, BEC in the sense of macroscopic occu- 
pation of the ground state, occurs when the short-range 
interparticle interactions are not too strong. Our results 
are in qualitative agreement with the recent PIMC sim- 
ulations |lj|] of similar systems. Both 2D and ID sim- 
ulation results indicate a phase transition for a 
finite number system, in contrast to the situation in the 
thermodynamic limit. Since systems of much larger size 
can be studied within the present approach, our work 
complements the PIMC calculations. 

The possibility of studying the tunneling phenomenon 
of condensed bosons in spatially different regions sepa- 
rated by a barrier has recently attracted some attention 
]3l| |34| . In particular, Dalfovo et al. |?2| have shown that 
a Josephson-type tunneling current may exist for bosons 
under the influence of a double-well trap potential. Za- 
pata et al. have estimated the Josephson coupling 
energy in terms of the condensate density. It is inter- 
esting to speculate on such a possibility in the present 
case, since the effective potential in our description is of 
the form of a double- well potential (cf. Fig. 2). In our 
treatment, the interaction effects modify the single-well 
trap potential into one which exhibits two minima. Thus 
if we think of this effective potential as the one seen by 
the condensed bosons and according to the general ar- 
guments |3l] |34|] based on two weakly connected systems 
we should have an oscillating flux of particles when the 
chemical potential in the two wells is different. Any con- 
figuration with Nl 7^ Nr which is always the case for 
odd number of bosons will result in an oscillatory mo- 
tion. It would be interesting to explore these ideas in 
future work. 
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FIG. 8. The temperature dependence of the specific heat 
Cv for various interaction strengths r\ and TV = 10 5 particles. 



IV. CONCLUDING REMARKS 

In this work we have applied the mean-field, semi- 
classical two-fluid model to interacting bosons in ID 
power-law trap potentials. We have found that for a 
range of interaction strengths the behavior of the thermo- 
dynamic quantities resembles to that of non-interacting 
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